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The  authors  propose  a  Ge/Ge0  76Si0  igSn0  05  quantum  cascade  laser  using  intersubband  transitions  at 
L  valleys  of  the  conduction  band  which  has  a  “clean”  offset  of  150  meV  situated  below  other  energy 
valleys  (r,X).  The  entire  structure  is  strain-free  because  the  lattice-matched  Ge  and 
Ge076Si0  19Sn0  05  layers  are  to  be  grown  on  a  relaxed  Ge  buffer  layer  on  a  Si  substrate.  Longer 
lifetimes  due  to  the  weaker  scattering  of  nonpolar  optical  phonons  reduce  the  threshold  current  and 
potentially  lead  to  room  temperature  operation.  ©  2007  American  Institute  of  Physics. 
[DOI:  10.1063/1.2749844] 


Electrically  pumped  Si-based  lasers  have  long  been 
sought  because  they  serve  as  light  sources  for  monolithic 
integration  of  Si  electronics  with  photonic  components  on 
the  same  Si  wafer.  Unfortunately,  Si  has  not  been  a  material 
of  choice  for  luminescence  applications  owing  to  its  indirect 
band  gap.  It  has  been  proposed  that  lasers  based  on  intersub¬ 
band  transitions  (ISTs)  in  SiGe  quantum  wells  (QWs)  could 
circumvent  the  issue  of  band  gap  indirectness.1  In  addition, 
SiGe  QWs,  being  a  nonpolar  material,  are  expected  to  have 
longer  intersubband  lifetimes  that  reduce  the  threshold  cur¬ 
rent  and  to  be  free  from  the  reststrahl  absorption  that  is  found 
in  III-V  quantum  cascade  lasers  (QCLs).  Various  groups 
have  obtained  electroluminescence  from  Si-rich  Si/SiGe 
quantum  cascade  structures, 2-4  but  lasing  has  eluded  re¬ 
searchers  up  to  now.  Those  efforts  all  have  one  scheme  in 
common:  holes  instead  of  electrons  are  used  for  1ST  because 
most  of  the  band  offset  between  the  Si-rich  SiGe  QWs  and 
the  Si  barriers  is  in  the  valence  band.  There  are  a  number  of 
difficulties  associated  with  valence  band  Si/SiGe  QCLs. 
First,  the  strong  mixing  of  heavy-hole,  light-hole,  and  split- 
off  bands  makes  the  QCL  design  exceedingly  cumbersome 
and  adds  a  great  degree  of  uncertainty.  Second,  the  large 
effective  mass  of  heavy  holes  hinders  carrier  injection  effi¬ 
ciency  and  leads  to  small  1ST  oscillator  strength  between 
laser  states.  Third,  for  any  significant  band  offset  needed  to 
implement  QCLs,  lattice-mismatch-induced  strain  in  SiGe 
QWs  is  likely  to  limit  total  structural  thickness  in  order  to 
avoid  generation  of  structural  defects.  Recently,  a  conduction 
intersubband  approach  was  proposed  to  construct  a  Ge/SiGe 
QCL  using  strained  Ge  QWs  and  SiGe  alloy  barriers.5  That 
structure  effectively  avoids  the  valence  band  complexity,  but 


“Electronic  mail:  hhcheng@ntu.edu. tw 


the  two  A2  valleys  along  the  (001)  growth  direction  are  still 
entangled  with  the  L  valleys  in  the  conduction  band,  leading 
to  design  complexity  and  potentially  creating  additional  non- 
radiative  decay  channels  for  the  upper  laser  state. 

In  this  letter,  we  propose  to  employ  Ge/Ge1_v_vSi(Snv 
heterostructures  to  develop  L-valley  QCLs.  Ge|_x_ySirSny  al¬ 
loys  have  been  studied  for  the  possibility  of  forming  direct 
band  gap  semiconductors.6  9  Since  the  initial  growth  of  this 
alloy,10  device-quality  epilayers  with  a  wide  range  of  alloy 
contents  have  been  achieved.  Incorporation  of  Sn  provides 
the  opportunity  to  engineer  separately  the  strain  and  band 
structure  since  we  can  vary  the  Si  (x)  and  Sn  (v)  composi¬ 
tions  independently.  Certain  alloy  compositions  of  this  mate¬ 
rial  system  offer  three  advantages:  (1)  the  possibility  of  a 
“cleaner”  conduction  band  lineup  in  which  the  L  valleys  in 
both  well  and  barrier  sit  below  other  valleys  (T.A-),  (2)  an 
electron  effective  mass  along  the  (001)  growth  direction  that 
is  much  lower  than  the  prior  heavy-hole  mass,  and  (3)  a 
lattice-matched  structure  that  is  entirely  strain-free.  In  addi¬ 
tion,  recent  advances  in  the  direct  growth  of  Ge  layer  on  Si 
provide  a  relaxed  matching  buffer  layer  on  a  Si  substrate 
upon  which  the  strain-free  QCL  is  grown,11  although  precise 
control  of  the  layer  thickness  and  interface  sharpness  is 
yet  to  be  accomplished  before  the  proposed  QCL  can  be 
implemented. 

Since  band  offsets  between  ternary  Sn-containing  alloys 
and  Si  or  Ge  are  not  known  experimentally,  we  have  calcu¬ 
lated  the  conduction  band  minima  for  a  lattice-matched  het¬ 
erostructure  consisting  of  Ge  and  a  ternary  Ge!_;c_>,SiASnv 
based  on  Jaros’  band  offset  theory,1”  which  is  in  good  agree¬ 
ment  with  experiment  for  many  heterojunction  systems.  For 
example,  this  theory  predicts  an  average  valence  band  offset, 
ALu  av=0.48  eV  for  a  Ge/Si  heterostructure  (higher  energy 
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FIG.  1.  Conduction  band  minima  at  L,T,X  points  of  Go,  ,  ,Si,Sn,  that  is 
lattice  matched  to  Ge. 

on  the  Ge  side),  close  to  the  accepted  value  of  A Evm 
=  0.5  eV.  The  basic  ingredients  of  our  calculation  are  the 
average  (between  heavy,  light,  and  split-off  hole  bands)  va¬ 
lence  band  offset  between  the  two  materials  and  the  compo¬ 
sitional  dependence  of  the  band  structure  of  the  ternary  alloy. 
For  a  Ge/a-Sn  interface  Jaros’  theory  predicts  A  £  av 
=  0.69  eV  (higher  energy  on  the  Sn  side).  For  the 
Ge  |  _x_,.SixSn,,/Ge  interface  we  have  used  the  customary  ap¬ 
proach  for  alloy  semiconductors,  interpolating  the  average 
valence  band  offsets  for  the  elementary  heterojunctions 
Ge/Si  and  Ge/a-Sn.  Thus  we  used  A  Evay(x,y) 
=£„jav(GeSiSn)-£i;  av(Ge)  =  -0.48x+0.69y  (in  eV).  Once  the 
average  valence  band  offset  is  determined,  the  energies  of 
individual  conduction  band  edges  in  the  Gei_l_vSi1.Sn>,  alloy 
can  be  calculated  relative  to  those  in  Ge  from  the  composi¬ 
tional  dependence  of  the  spin-orbit  splitting  of  the  top  va¬ 
lence  band  states  and  the  compositional  dependence  of  the 
energy  separations  between  those  conduction  band  edges  and 
the  top  of  the  valence  band  in  the  alloy.  '  We  have  assumed 
that  all  required  alloy  energies  can  be  interpolated  between 
the  known  values  for  Si,  Ge,  and  a-Sn  as 

EGcSiSn(x,y)  =  EGe(l  -x-  y)  +  Esix 

+  ESny  ~  bGeSi(\  -x-y)x 

-  bGcSn(  1  -  x  -  >')}’  -  bSiSltxy .  ( 1 ) 

The  bowing  parameters  bGe si,  bGeSn,  and  Z?siSn  have  been  dis¬ 
cussed  in  Refs.  14  and  15.  Finally,  for  the  indirect  conduc¬ 
tion  band  minimum  near  the  X  point,  Weber  and  Alonso  find 
£x=0.931  +  0.018.r+0.206.r2  (in  eV)  for  Ge,_vSiv  alloys.16 
On  the  other  hand,  the  empirical  pseudopotential  calculations 
of  Cohen  and  Bergstresser  place  this  minimum  at  0.90  eV  in 
cr-Sn,  virtually  the  same  as  its  value  in  pure  Ge.17  We  thus 
assume  that  the  position  of  this  minimum  in  ternary 
Ge^^SijSm  alloys  is  independent  of  the  Sn  concentration 
y.  This  is  qualitatively  in  agreement  with  the  virtual-crystal 
calculations  of  Jenkins  and  Dow,7  which  find  that  the 
A- related  minima  have  the  weakest  compositional  depen¬ 
dence.  The  conduction  band  minima  results  are  shown  in 
Fig.  1  for  Sn  concentrations  0<y<  0.1.  The  Si  concentration 
x  was  calculated  using  Vegard’s  law  in  such  a  way  that  the 
ternary  Ge1_r_vSi,.Snv  is  exactly  lattice  matched  with  Ge. 

It  can  be  seen  from  Fig.  1  that  a  conduction  band  offset 
of  150  meV  at  L  valleys  can  be  obtained  between  lattice- 
matched  Ge  and  Ge0  76Si0 ,i9Sn0  05  alloy  while  all  other  con¬ 
duction  band  valleys  (r,A,  etc)  are  above  the  L- valley  band 
Downloaded  30  Jul  2007  to  146.153.144.35.  Redistribution  subject 
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FIG.  2.  (Color  online)  L-valley  conduction  band  profile  and  squared  enve¬ 
lope  functions  under  an  electric  field  of  lOkV/cm.  Layer  thicknesses  in 
angstrom  are  marked  with  bold  numbers  for  Ge  QWs  and  regular  for  Ge- 
SiSn  barriers.  Array  marks  the  injection  barrier. 

edge  of  the  Ge0  76Si0  19Sn0  05  barrier.  This  band  alignment 
presents  a  desirable  alloy  composition  from  which  a  QCL 
operating  at  L  valleys  can  be  designed  using  Ge  as  QWs  and 
Ge076Sio.i9Sn0.o5  as  barriers  without  the  complexity  arising 
from  other  energy  valleys. 

Figue  2  shows  the  QCL  structure  based  upon 
Ge/Ge076Si0  19Sn005  QWs.  Only  L-valley  conduction  band 
lineups  are  shown  in  the  potential  diagram  under  an  applied 
electric  field  of  10  kV/cm.  In  order  to  solve  the  Schrodinger 
equation  to  yield  subbands  and  their  associated  envelope 
functions,  it  is  necessary  to  determine  the  effective  mass  m, 
along  the  (001)  growth  direction  (z)  within  the  constant- 
energy  ellipsoids  at  the  L  valleys  along  the  (111)  direction 
which  is  tilted  with  respect  to  (001).  Using  the  L-valley  prin¬ 
cipal  transverse  effective  mass  in,  =  0.08mo  and  the  longitu¬ 
dinal  effective  mass  w/=L60mo,  we  obtain  m.={2l3mt 
+  l/3m/)_1=0.12mo  where  m0  is  the  free  electron  mass.  The 
squared  magnitudes  of  all  envelope  functions  are  plotted  at 
energy  positions  of  their  associated  subbands.  As  shown  in 
Fig.  2,  each  period  of  the  QCL  has  an  active  region  for  lasing 
emission  and  an  injector  region  for  carrier  transport.  These 
two  regions  are  separated  by  a  30  A  barrier.  The  active  re¬ 
gion  is  constructed  with  three  coupled  Ge  QWs  that  give  rise 
to  three  subbands  marked  1,  2,  and  3.  The  lasing  transition  at 
the  wavelength  of  49  /im  is  between  the  upper  laser  state  3 
and  the  lower  laser  state  2.  The  injector  region  consists  of 
four  Ge  QWs  of  decreasing  well  widths  all  separated  by 
20  A  Ge0  76Si0  i9Sn0  05  barriers.  The  depopulation  of  lower 
state  2  is  through  scattering  to  state  1  and  to  the  miniband 
downstream  formed  in  the  injector  region.  These  scattering 
processes  are  rather  fast  because  of  the  strong  overlap  be¬ 
tween  the  involved  states.  Another  miniband  in  the  injector 
region  formed  of  quasibound  states  is  situated  45  meV  above 
the  upper  laser  state  3,  effectively  preventing  escape  of  elec¬ 
trons  from  upper  laser  state  3  into  the  injector  region. 

The  nonradiative  transition  rates  between  different  sub¬ 
bands  in  such  a  low-doped  nonpolar  material  system  with 
low  injection  current  should  be  dominated  by  deformation- 
potential  scattering  of  nonpolar  optical  and  acoustic 
phonons.  Because  of  the  large  energy  separation  between  L 
and  T  valleys,  and  between  L  and  X  valleys,  the  L-F  and  L-X 
intervalley  scattering  processes  are  negligible.  Between  the 
four  equivalent  L  valleys,  it  can  be  argued  through  the  matrix 

element  for  electron-phonon  interaction  that  the  inter 
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FIG.  3.  Lifetimes  and  threshold  current  density  as  a  function  of  operating 
temperature. 


L- valley  scattering,  an  umklapp  process,  is  much  weaker 
than  the  intra-L-valley  scattering.  The  lifetimes  are  therefore 
determined  by  the  intra-L-valley  scattering.  For  this  Ge-rich 
structure,  we  have  used  bulk-Ge  phonons  for  calculation  of 
the  scattering  rate  to  yield  lifetimes  for  the  upper  laser  state 
t3  and  the  lower  laser  state  t2,  as  well  as  the  3— >2  scattering 
time  t32.18  They  are  shown  in  Fig.  3  as  a  function  of  operat¬ 
ing  temperature.  These  lifetimes  are  at  least  one  order  of 
magnitude  longer  than  those  of  III-V  QCLs  owing  to  the 
nonpolar  nature  of  GeSiSn  alloys.  The  necessary  condition 
for  population  inversion  r32  >  r2  is  satisfied  throughout  the 
temperature  range.  Using  these  predetermined  lifetimes  in 
the  population  rate  equation  under  current  injection: 


dN3  7]J  N3  -  N-$ 

dt  e  r3 


dNs 

dt 

V. 


Ns- Ns 

t32 


n2-n2 


(2) 


where  e  is  the  electron  charge,  Ni(i=2,3)  is  the  area  carrier 
density  per  period  in  subband  i  under  injected  current  density 
J  with  an  injection  efficiency  77,  and  Nt  is  the  area  carrier 
density  per  period  due  to  thermal  excitation  of  n  doping. 
Solving  the  above  rate  equation  at  steady  state  yields  popu¬ 
lation  inversion 


Ns-N2=Ts\  1 


^32/  e 


■(N2~Ns), 


(3) 


which  can  be  then  used  to  evaluate  the  optical  gain  of  the 
TM-polarized  mode  as19 


2e2A\(3\p,\2)\2 

s0cnm2.yLp{ha)L) 


—)--(n2-nj 

t32/  e 


(4) 


where  ea  is  the  permittivity  in  vacuum,  c  the  speed  of  light  in 
vacuum,  h  the  Planck  constant,  and  (3[p.|2)  the  momentum 
matrix  between  the  two  laser  states.  Other  parameters  are  as 
follows:  index  of  refraction  /?  =  3.97,  lasing  transition  energy 
ha>L=25  meV,  full  width  at  half  maximum  y=  1 0  meV, 
length  of  one  period  of  the  QCL  Lp= 532  A,  area  doping 
density  per  period  of  1010/cm2,  and  unit  injection  efficiency 
V=  1  ■ 


We  have  simulated  the  TM-polarized  mode  in  a  QCL 
structure  of  40  periods  that  is  confined  by  double-Au- 
plasmon  waveguide  and  obtained  near  unity  optical  confine¬ 
ment  T~  1.0  and  waveguide  loss  aw=  1 10/cm.  Assuming  a 
mirror  loss  am=  10/cm  for  a  typical  cavity  length  of  1  mm, 
the  threshold  current  density  /th  can  be  calculated  from  the 
balancing  relationship,  T gth=  aw+  am.  The  result  is  shown  in 
Fig.  3  for  /th  that  ranges  from  22  A/cm2  at 
5  K  to  550  A/cm2  at  300  K.  These  threshold  values  are 
lower  than  those  of  III-V  QCLs  as  a  result  of  the  longer 
scattering  times  due  to  nonpolar  optical  phonons. 

In  summary,  we  propose  a  Ge/Ge076Si0.i9Sn0  05  QCL 
that  operates  at  L  valleys  of  the  conduction  band.  According 
to  our  estimation  of  the  band  lineup,  this  particular  alloy 
composition  gives  a  clean  conduction  band  offset  of 
150  meV  at  L  valleys  with  all  other  energy  valleys  conve¬ 
niently  out  of  the  way.  All  QCL  layers  are  lattice  matched  to 
a  Ge  buffer  layer  on  a  Si  substrate  and  the  entire  structure  is 
therefore  strain-free.  The  electron  effective  mass  along  the 
growth  direction  is  much  lighter  than  that  of  heavy  holes 
bringing  a  significant  improvement  in  tunneling  rates  and 
oscillator  strengths.  The  lasing  wavelength  of  this  device  is 
49  yum.  With  different  GeSiSn  alloy  compositions  that  are 
lattice  matched  to  Ge,  QCLs  can  be  tuned  to  lase  at  other 
desired  wavelengths.  Lifetimes  determined  from  the 
deformation-potential  scattering  of  nonpolar  optical  and 
acoustic  phonons  are  at  least  an  order  of  magnitude  longer 
than  those  in  III-V  QCLs  with  polar  optical  phonons,  leading 
to  a  reduction  in  threshold  current  density  and  the  possibility 
of  room  temperature  operation. 

This  work  is  supported  in  part  by  the  Air  Force  Office  of 
Scientific  Research. 
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